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Combined constraints from the CDF and DO Collaborations on models of the Higgs boson with 
exotic spin J and parity P are presented and compared with results obtained assuming the standard 
model value J F = 0 + . Both collaborations analyzed approximately 10 fb” 1 of proton-antiproton 
collisions with a center-of-mass energy of 1.96 TeV collected at the Fermilab Tevatron. Two models 
predicting exotic Higgs bosons with J p = 0 _ and J P = 2 + are tested. The kinematic properties of 
exotic Higgs boson production in association with a vector boson differ from those predicted for the 
standard model Higgs boson. Upper limits at the 95% credibility level on the production rates of 
the exotic Higgs bosons, expressed as fractions of the standard model Higgs boson production rate, 
are set at 0.36 for both the J p = 0“ hypothesis and the J p = 2 + hypothesis. If the production 
rate times the branching ratio to a bottom-antibottom pair is the same as that predicted for the 
standard model Higgs boson, then the exotic bosons are excluded with significances of 5.0 standard 
deviations and 4.9 standard deviations for the J p = 0“ and J p = 2 + hypotheses, respectively. 

PACS numbers: 13.85.Rm, 14.80.Bn, 14.80.Ec 


The Higgs boson discovered by the ATLAS [l]] and 
CMS Q Collaborations in 2012 using data produced in 
proton-proton collisions at the Large Hadron Collider 
(LHC) at CERN allows many stringent tests of the elec- 
troweak symmetry breaking in the standard model (SM) 
and extensions to the SM to be performed. To date, 
measurements of the Higgs boson’s mass and width [§- 
HI, its couplings to other particles B-flll . and its spin 
and parity quantum numbers J and P 10-|l6| are consis¬ 
tent with the expectations for the SM Higgs boson. The 
CDF and DO Collaborations at the Fermilab Tevatron 
observed a 3.0 standard deviation (s.d.) excess of events 


consistent with a Higgs boson signal, largely driven by 
those channels sensitive to the decay of the Higgs boson 
to bottom quarks (H —i bb) 00. The Tevatron data 
are also consistent with the predictions for the properties 
of the SM Higgs boson [Tsl - E^ . 

Ref. (23| proposed to use the Tevatron data to test 
models for the Higgs boson with exotic spin and parity, 
using events in which the exotic Higgs boson X is pro¬ 
duced in association with a W or a Z boson and decays to 
a bottom-antibottom quark pair, X —> bb. This proposal 
used two of the spin and parity models in Ref. [24| . one 
with a pseudoscalar J p = 0 _ state and the other with 
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a graviton-like J p = 2 + state. For the SM Higgs bo¬ 
son, which has J p = 0 + , the differential production rate 
near threshold is linear in /3, where f) = 2p/y/l, p is the 
momentum of the X boson in the VX (V = W or Z) ref¬ 
erence frame, and y/§ is the total energy of the VX system 
in its rest frame. For the pseudoscalar model, the depen¬ 
dence is proportional to /3 3 . For the graviton-like model, 
the dependence is proportional to /3 5 - however, not all 
J p = 2 + models share this /3 5 factor .23]. These powers 
of /3 alter the kinematic distributions of the observable 
decay products of the vector boson and the Higgs-like 
boson X, most notably the invariant mass of the VX 
system, which has a higher average value in the J p = 0“ 
hypothesis than in the SM 0 + case, and higher still in 
the J p = 2 + hypothesis. These models predict neither 
the production rates nor the decay branching fractions of 
the A' particles. 

The ATLAS and CMS Collaborations recently re¬ 
ported strong evidence for Higgs boson decays to 
fermions [25 |t 30| , with sensitivity dominated by the H —> 
t + t~ decay mode, though they have not yet performed 
spin and parity tests using fermionic decays. The par¬ 
ticle decaying fermionically for which the Tevatron also 
found evidence might not be the same as the particle dis¬ 
covered through its bosonic decays at the LHC. Tests of 
the spin and parity frl with Tevatron data therefore 
provide unique information on the identity and prop¬ 
erties of the new particle or particles. The CDF and 
DO Collaborations have re-optimized their SM Higgs bo¬ 
son searches to test the exotic Higgs boson models in 


the WH 


£vbb 


ZH 


£ + £ bb 33j, |34|, and 


WH + ZH —> $ T bb [35], ]36( channels, where £ = e or p 
and lft T is the missing transverse energy 3J. In this let¬ 
ter we report a combination of the CDF [21] and DO [22] 


studies of the J p assignments of the state X, with mass 
mx = 125 GeV/c 2 , in the X —y bb decay. 

The CDF and DO detectors are multipurpose solenoidal 
spectrometers surrounded by hermetic calorimeters and 
muon detectors designed to study the products of 1.96 
TeY proton-antiproton (pp) collisions 


|38 


)• All 

searches combined here use the complete Tevatron data 
sample, which, after data quality requirements, corre¬ 
sponds to 9.45 - 9.7 fb -1 of integrated luminosity, de¬ 
pending on the experiment and the search channel. 

Standard model Higgs boson signal events are sim¬ 
ulated u sing the leading-order (LO) calculation from 
pythia [40], with CTEQ5L (CDF) and CTEQ6L1 
(DO) f4l| parton distribution functions (PDFs). The 
J p = 0” and J p = 2 + signal samples are generated 
using MADGRAPH 5 version 1.4.8.4 42]^ with modifica¬ 
tions provided by the authors of Ref. [23]. Subsequent 
particle showering is modeled by PYTHIA. We normalize 
the SM Higgs boson rate predictions to the highest-order 
calculations available. The WH and ZH cross sections are 
calculated at next-to-next-to-leading-order (NNLO) pre¬ 
cision in the strong interaction, and next-to-leading-order 


(NLO) precision in the electroweak corrections I l3l-| l6l. 
We use the branching fractions for Higgs boson decay 
from Ref. [47}. These rely on calculations using hde- 
cay 0 and prophecy4f fUj. 

The predictions of the dominant background rates and 
kinematic distributions are treated in the following way. 
Diboson ( WW, WZ , and ZZ) Monte Carlo (MC) sam¬ 
ples are normalized using the NLO calculations from 
MCFM j50]. For ti. we use a production cross section of 
7.04 ± 0.70 pb 51], which is based on a top quark mass 
of 173 GeV/c 2 [H3 and MSTW 2008 NNLO PDFs [H3]. 
The single top quark production cross section is assumed 
to be 3.15 ± 0.31 pb (54]. For details of the generators 
used, see Ref. [55] . Data-driven methods are used to nor¬ 
malize the V plus light-flavor and heavy-flavor jet back¬ 
grounds [60] using V data events containing no b- ta gge d 
jets [glj, which have negligible signal content 62, |63( . 
The MC modeling of the kinematic distributions of the 
background predictions is described in Refs. [3ll436j . 

The event selections are similar (CDF), or identical 
(DO), to those used in their SM counterparts 31-361]. 
For the WH —> £vbb analyses, events are selected with 
one identified lepton (e or p), jets, and large $ T . For 
the CDF WH —> £vbb analysis, only events with two jets 
are used. Events are classified into separate categories 
based on the quality of the identified lepton. Separate 
categories are used for events with a high-quality muon 
or central electron candidate, an isolated track, or a for¬ 
ward electron candidate. Within the lepton categories, 
five exclusive 6-tag categories, comprising two single-tag 
and three double-tag categories, are formed. The mul¬ 
tivariate 6-tagger used by CDF 64| was trained on SM 
Higgs boson signal MC events. Few of these events con¬ 
tained jets with with transverse energy Ex > 200 GeV 
and thus the tagger does not perform well for such jets. 
Hence, only jets with Et < 200 GeV are considered. For 
the DO WH —> £vbb analysis, events are selected with 
two or three jets. The data are split by lepton flavor 
and jet multiplicity (two or three jet subchannels), and 
by the output of the 6-tagging algorithm applied to all 
selected jets in the event. This channel, along with the 
other two DO channels, uses a multivariate 6-tagging al¬ 
gorithm 65, [66|. Four exclusive 6-tag categories, one 
single-tag and three double-tag, are formed. In the SM 
Higgs boson search, boosted decision trees are used as 
the final discriminating variables; here they are used to 
further subdivide the selected data sample into high- and 
low-purity categories. 

The ZH —> £ + £~bb analyses require two isolated lep¬ 
tons and at least two jets. The CDF analysis separates 
events into one single- and three double-6-tag samples 
and uses neural networks to select loose dielectron and 
dimuon candidates. The jet energies are corrected for 
$ T using a neural network ffjjj]- The CDF analysis uses 
a multistep discriminant based on neural networks, where 
two discriminant functions are used to define three sep- 
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arate regions of the final discriminant function. The 
DO ZH —> £ + £~bb analysis separates events into non¬ 
overlapping samples of events with either a single or dou¬ 
ble b- tag. To increase the signal acceptance, the selection 
criteria for one of the leptons are loosened to include iso¬ 
lated tracks not reconstructed in the muon detector and 
electron candidates from the intercryostat region of the 
DO detector. Combined with the dielectron and dimuon 
categories, these provide four independent lepton sub¬ 
channels. A kinematic fit is used to optimize reconstruc¬ 
tion. Random forests (RF) of decision trees EH EH are 
used to provide the final variables in the SM Higgs boson 
search. The first RF is designed to discriminate against tt 
events and divides events into tt-enriched and tt-depleted 
single-tag and double-tag regions. Only events in the 
tt-depleted regions are considered in this study. These 
regions contain approximately 94% of the SM signal. 

For the ZH —> vvbb analyses, the selections used by 
CDF and DO are similar to the WH selections, except 
that all events with isolated leptons are rejected and more 
stringent techniques are applied to reject the multijet 
background. In a sizable fraction of WH —> £vbb sig¬ 
nal events, the lepton is undetected. Such events often 
are selected in the ZH — > vvbb samples, so these analyses 
are also referred to as VH — > $ T bb. The CDF analysis 
uses three non-overlapping 6-tag categories (two double- 
and one single-tag) and two jet categories (two- or three- 
jet events), giving a total of six subchannels. In the DO 
analysis, exactly two jets are required and two exclusive 
double-tag categories are defined using the sum of the 
6-tagging outputs for each of the two selected jets. 

Both CDF and DO have a 50% larger acceptance for the 
J p = 0“ and J p = 2 + signals in the ZH — > vvbb analyses 
compared with the SM Higgs boson signal, largely due 
to the fact that the exotic signal events are more likely 
to pass the trigger thresholds for $ T , a consequence of 
the larger average VX invariant masses. The other two 
channels, WH —t £vbb and ZH —> £ + t~bb, do not benefit 
as much from the additional $ T in these events, as they 
rely on the lepton triggers, which are more efficient than 
the triggers in the relevant kinematic regions. 

Unlike their SM counterparts, these analyses are op¬ 
timized to distinguish the J p = 0~ and the J p = 2 + 
hypotheses from the SM 0 + hypothesis. The exotic par¬ 
ticles are considered either in addition to, or replacing, 
the SM Higgs boson. A mixture of all three states is not 
considered. 

The CDF multivariate analysis (MVA) discriminants 
were newly trained to separate the exotic Higgs boson 
signals from the SM backgrounds [21] . In the WH —> 
tvbb and VH —> $ T bb channels, events classified as 
background-like by the new discriminants are then clas¬ 
sified according to the SM-optimized MVA discriminants 
in order to improve the performance of tests between the 
SM and exotic hypotheses. 

Depending on the channel, DO uses either the recon¬ 


structed dijet mass or the MVA used in the SM Higgs bo¬ 
son search to separate events into high- and low-purity 
samples. The mass of the VX system is then used to 
discriminate between the exotic and SM hypotheses 221. 
For the ZH —> Ubb analysis the invariant mass of the two 
leptons and the two highest pr jets is used. For the ivbb 
and vvbb final states the transverse mass Mp is used, 
where Mj. = {Ej, + Ep ) 2 — (pp +Pr) 2 an d the trans¬ 
verse momenta of the Z and W bosons are taken to be 
Pp = |? T and pft = $ T + Pp, respectively. 


The number of contributing channels is large, and their 
sensitivities vary from one to another. To visualize the 
data in a way that emphasizes the sensitivity to the ex¬ 
otic signals, we follow Ref. [Til ]. Bins of the final discrim¬ 
inant for all channels are ordered by increasing signal- 
to-background ratio (s/6) and are shown in comparison 
with predicted yields from signal and background pro¬ 
cesses for the J p = 0“ and J p = 2 + searches in Fig. U| 
separately. The backgrounds are fit to the data in each 
case, allowing the systematic uncertainties to vary within 
their a priori constraints. The exotic signals are normal¬ 
ized to the SM cross section times branching ratio mul¬ 
tiplied by an exotic-signal scaling factor, p e xotic- They 
are shown in Fig. Q] with /i eX otic = 1- The scaling fac¬ 
tor for the SM Higgs boson signal is denoted by psm- A 
value of one for either psm or p e xotic corresponds to a 
cross section times branching ratio as predicted for the 
SM Higgs boson. Both figures show agreement between 
the background predictions and the observed data over 
five orders of magnitude with no evidence for an excess 
of exotic signal-like candidates. 

We follow Ref. fl8j and perform both Bayesian and 
modified frequentist calculations of the upper limits on 
exotic X boson production with and without SM Higgs 
production, best-fit cross sections allowing for the simul¬ 
taneous presence of a SM Higgs boson and an exotic X 
boson, and hypothesis tests for signals assuming various 
production rate times branching ratio values for the ex¬ 
otic bosons. Both methods use likelihood calculations 
based on Poisson probabilities that include SM back¬ 
ground processes and signal predictions for the SM Higgs 
and exotic bosons multiplied by their respective scaling 
factors, psm and /r eX otic- Systematic uncertainties on the 
predicted rates and on the shapes of the distributions and 
their correlations are treated as described in Ref. jl8j |. 
Theoretical uncertainties in cross sections and branching 
ratios are considered fully correlated between CDF and 
DO, and between analysis samples. The uncertainties on 
the measurements of the integrated luminosities, which 
are used to normalize the expected signal yields and the 
MC-based background rates, are 6.0% (CDF) and 6.1% 
(DO). Of these values, 4% arises from the inelastic pp 
cross section [70], which is fully correlated between CDF 
and DO. The dominant uncertainties on the backgrounds 
are constrained by the data in low s/6 regions of the 
discriminant distributions. Different methods were used 






FIG. 1: (color online). Distribution of log 10 (s/6) for CDF 
and DO data from all contributing search channels, for (a) the 
J p = 0 _ search and (b) the J p = 2 + search. The data are 
shown with points, and the expected exotic signals are shown 
with /iexotic = 1 stacked on top of the fitted backgrounds. The 
solid lines denote the predictions for the SM Higgs boson, and 
are not stacked. Underflows and overflows are collected into 
the leftmost and rightmost bins, respectively. 


by CDF and DO to estimate U+jets and multijet back¬ 
grounds and so their uncertainties are considered uncor¬ 
related. Similarly, the uncertainties on the data-driven 
estimates of the 6-tag efficiencies are considered uncorre¬ 
lated between CDF and DO, as are the uncertainties on 
the jet energy scales, the trigger efficiencies, and lepton 
identification efficiencies. We quote Bayesian upper lim¬ 
its and best-fit cross sections assuming uniform priors for 
non-negative signal cross sections, and we use the modi¬ 
fied frequentist method to perform the hypothesis tests. 
Systematic uncertainties are parameterized by nuisance 
parameters with Gaussian priors, truncated so that no 
predicted yield for any process in any search channel is 


negative. 

For both the J p = 0“ and J p = 2 + models, we com¬ 
pute two 95% credibility upper limits on /j ex ot i c , one as¬ 
suming /jsm = 1 and the other assuming /jsm = 0. The 
expected limits are the median expectations assuming 
no exotic boson is present. The results are listed in Ta¬ 
ble HI Two-dimensional credibility regions, which are the 
smallest regions containing 68% and 95% of the poste¬ 
rior probabilities, are shown in Fig. [2] The points in the 
(msm, Mexotic) planes that maximize the posterior prob¬ 
ability densities are shown as the best-fit values. These 
best-fit values are (/xsm=L 0, p 0 - = 0) for the search 
for the J p = 0“ state, and (/xsm=L1, P- 2 + = 0) for the 
search for the J p = 2 + state. We also derive upper limits 
on the fraction f JP = Atexotic/tpexotic+AGM), as functions 
of the total p = /x eX otic + A*sm, assuming a uniform prior 
probability density in non-negative fjp, extended to in¬ 
clude fractions larger than 1.0 in order not to saturate 
the limits at fjp = 0.95 for p < 0.6, where the test is 
weak. The results are shown in Fig. [3] 

In the modified frequentist approach 00 we com¬ 
pute p values for the discrete two-hypothesis tests, 
the SM Higgs boson hypothesis (the “null” hypothe¬ 
sis) (/xsm=1, l*exotic=0) and the exotic (“test”) hypoth¬ 
esis (/xsm— 0, /r eX otic=l), both assuming that SM back¬ 
ground processes are present. The choice of setting 
Mexotic = 1 in the test hypothesis is arbitrary; the sen¬ 
sitivity of the test is reduced if a smaller value is as¬ 
sumed. We use the log-likelihood ratio, LLR, defined 
to be —2 ln(p(data|test)/p(data|null), where the numer¬ 
ator and denominator are maximized over systematic 
uncertainty variations 0. The LLR distributions are 
shown in the supplemental material [73]. We define 
the p values p nu n = P(LLR < LLR 0 b s |SM) and p tes t = 
P(LLR > LLR 0 b s (exotic). The median expected p values 
Tnunimed in the test hypothesis andp™ med in the SM hy¬ 
pothesis quantify the sensitivities of the two-liypothesis 
tests for exclusion and discovery, respectively. Table [TT] 
lists these p values for both exotic models, as well as 
CL S = ptest/(I — Pnuii) 0] for the Tevatron combina¬ 
tion. To compute pt es t and the expected values of p nu ii 
and ptest, Wilks’s theorem is used [741 ] . 

The similarity of the limits and p values obtained for 
the J p = 0“ and the J p = 2 + searches is expected since 
the exotic models predict excesses in similar portions of 
kinematic space. 

In summary, we combine CDF’s and DO’s tests for the 
presence of a pseudoscalar Higgs boson with J p = 0~ 
and a graviton-like boson with J p = 2 + in the WX —y 
Ivbb, the ZX —> £ + £~bb, and the VX —> lf T bb search 
channels using models described in Ref. 23]. The masses 
of the exotic bosons are assumed to be 125 GeV/c 2 . No 
evidence is seen for either exotic particle, either in place 
of the SM Higgs boson or produced in a mixture with a 
J p = 0+ Higgs boson. In both searches, the best-fit cross 
section times the decay branching ratio into a bottom- 
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TABLE I: Observed and median expected Bayesian upper 
limits at the 95% credibility level on /.i eX otic for the pseu¬ 
doscalar ( J p = 0 _ ) and graviton-like (J p = 2 + ) boson mod¬ 
els, assuming either that the SM Higgs boson is also present 
(/tsm = 1) or absent (/xsm = 0). 


Tevatron Run II, L < 10 fb 




o 

zL 


0.8 


Channel 


Observed 

(Limit/o- SM ) 

Median Expected 
(Limit/crsM) 

0.6 


' 

J P = 0”, /J.SM 

= 0 

0.36 

0.32 

0.4 


_ 

J P = 0 _ , pLSM 

= 1 

0.29 

0.32 



- 

J P = 2 + , /iSM 

= 0 

0.36 

0.33 




J P = 2 + , /iSM 

= 1 

0.31 

0.34 

0.2 




“T 


VX^ Vbb 


• Best fit 
aSM 

□ 68% C.L. jP =0 - (a) 

I 95% C.L. 


TABLE II: Observed (obs) and median expected (med) LLR. 
values and p values for the combined CDF and DO searches 
for the pseudoscalar ( J p = 0“) boson and the graviton-like 
(J p = 2 + ) boson. The p values are listed, and the corre¬ 
sponding significances in units of standard deviations, using 
a one-sided Gaussian tail calculation, are given in parenthe¬ 
ses. The two hypotheses tested are (/tsm, Mexotic) = (1,0) and 
(0,1) for the SM and the exotic models, respectively. 


Analysis 

Hi 

II 

o 

1 

J p = 2+ 

LLRobs 

27.1 

25.7 

T T pSM 

LLn med 

23.7 

21.8 

t t t) exotic 

J-jTjiUmed 

-29.9 

-29.6 

Pnull 

0.63 (-0.34) 

0.66 (-0.41) 

exotic 
/''null,med 

1.8xl0 -8 (5.5) 

1.9xl0" 8 (5.5) 

Ptest 

9.4xl0 -8 (5.2) 

1.9xl0 -7 (5.1) 

Ptert, med 

4.7xl0 -7 (4.9) 

1.2xl0 -6 (4.7) 

CL S 

2.6xl0 -7 (5.0) 

5.6xl0 -7 (4.9) 

/''it SM 
'^■^ J s,med 

9.4xl0 -7 (4.8) 

2.3xl0 -6 (4.6) 


1-1 


SM 


CM 


1 

0.8 

0.6 

0.4 

0.2 

0 


Tevatron Run II, L < 10 fb 


• BesTfit 
aSM 

□ 68% C.L. 
B95% C.L. 


0 


J p = 2" 1 


VX^ Vbb 

(b) 


V, 


SM 


FIG. 2: (color online). Two-dimensional credibility regions 
in the (/iexotic, psm) plane, for the combined CDF and DO 
searches for (a) the pseudoscalar (J p = 0“) boson, and (b) 
the graviton-like (J p = 2 + ) boson. 


antibottom quark pair of a J p = 0 + signal component 
is consistent with the prediction of the SM Higgs boson. 
The Bayesian posterior probability densities for_the J p = 
CP and J p = 2 + searches are shown in Ref. [13]. 

Upper limits at 95% credibility on the rate of the pro¬ 
duction of an exotic Higgs boson in the absence of a SM 
J p = 0 + signal are set at 0.36 times the SM Higgs pro¬ 
duction rate for both the J p = 0“ and the J p = 2 + 
hypotheses. If the production rate of the hypothetical 
exotic particle times its branching ratio to a bottom- 
antibottom quark pair is the same as that predicted for 
the SM Higgs boson, then the exotic models are excluded 
with significances of 5.0 s.d. and 4.9 s.d. for the J p = 0” 
and J p = 2 + hypotheses, respectively. 
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FIG. 4: Posterior probability density distributions for the combined searches for exotic J p = 0~ and J p = 2 + bosons, (a) 
The posterior probability density as a function of /r 0 - assuming /ism = 1 and a uniform prior density for non-negative fj, 0 - , 
(b) the posterior probability density as a function of fj, 0 - assuming fi sm = 0 and a uniform density for non-negative /r 0 -, 
and (c) the posterior probability density as a function of the fraction of exotic boson production, + A 4 sm), assuming 

/r 0 - + Msm = 1, and a uniform prior density for non-negative values of the fraction. The dashed vertical lines indicate the 
observed limits. Figures (d), (e), and (f) show the corresponding results for the J p = 2 + boson searches. 
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FIG. 5: Distributions of LLR for the combined CDF and DO searches for (a) the pseudoscalar (J p = 0 _ ) boson, and (b) the 
graviton-like (J p = 2 + ) boson. The LLR distributions are shown separately assuming that an exotic particle is present with 
/^exotic = 1 plus SM backgrounds, and if the SM Fliggs boson plus SM backgrounds are present. The observed values of LLR 
are shown with vertical lines. Shaded regions show the 68% and 95% confidence level regions on the distributions assuming the 
SM Higgs boson is present, centered on the median expectation. 
































